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Abstract. Well-posedness for the initial value problem for a self-gravitating 
elastic body with free boundary in Newtonian gravity is proved. In the material 
frame, the Euler-Lagrange equation becomes, assuming suitable constitutive 
properties for the elastic material, a fully non-linear elliptic-hyperbolic system 
with boundary conditions of Neumann type. For systems of this type, the 
initial data must satisfy compatibility conditions in order to achieve regular 
solutions. Given a relaxed reference configuration and a sufficiently small 
Newton's constant, a neigborhood of initial data satisfying the compatibility 
conditions is constructed. 



1. Introduction 

In Newtonian physics, the two-body problem is solvable for point particles mov- 
ing around their common center of gravity on Kepler ellipses. However, if one 
considers extended bodies, the situations changes drastically. Assuming that a so- 
lution exists for extended bodies, one can show that the centers of mass of the 
bodies move as point particles, but existence could not be established for a long 
time. 

The first sucessful attack on the problem was made by Leon Lichtenstein [21] 
who considered self-gravitating fluid bodies moving on circles about their center of 
gravity. In this case the Euler equations for a self-gravitating system become time 
independent in a coordinate system co-moving with the bodies. For the case of 
small bodies or widely separated large bodies, Lichtensten showed the existence of 
such solutions. 

Well-posednedness for the Cauchy problem for fluid bodies with free boundary 
was proved only recently. Lindblad [22] proved well-posedness for a non-relativistic 
compressible liquid body (i.e. positive boundary density) with free boundary. In 
this paper one can also find references to earlier work. A different proof of the 
result of Lindblad that is valid for both the relativistic and non-relativistic cases 
was given by Trakhinin [33] . 

The more singular case of fluids with vanishing boundary density is discussed in 
[12, 17]. Unfortunately, the problem of proving well-posedness for self-gravitating 
compressible fluid bodies with free boundary is still open in general relativity and 
even in Newtonian gravity. However, see [23] for the case of self-gravitating in- 
compressible fluids. See also [29, 10, 18] for results dealing with various restricted 
versions of the Cauchy problem for fluid bodies in general relativity. 

One can argue that the Cauchy problem should be simpler to handle if the bodies 
consist of elastic material, since such a body can, at least in the case of small bodies, 
be said to 'have a shape of its own'. In this paper we solve the boundary initial value 
problem for self-gravitating deformations of a relaxed body. In particular, consider 
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a relaxed (i.e. with vanishing stress) elastic body in the absence of gravity. Then 
we show existence of a self-gravitating solution for initial data close to those of the 
static relaxed body and for a small gravitational constant. Thus, the motion of the 
body will consist of small nonlinear oscillations around the equilibrium solution. 

The Cauchy problem for a Newtonian elastic body with free boundary is, under 
suitable constitutive assumptions on the elastic material, a fully non-linear elliptic- 
hyperbolic initial-boundary value problem with boundary values of Neumann type. 
The case of purely hyperbolic problems of this type is covered by a theorem of Koch 
[19]. The method of proof of this theorem can be adapted to include the non-local 
elliptic terms in the which appear in the system of differential equations considered 
in this paper via the Newtonian gravitational field. 

In order to achieve a regular evolution for the initial-boundary value problem, 
the initial data must satisfy certain compatibility conditions. The problem of con- 
structing an open neighborhood of initial data satisiying the compatibility condi- 
tions given one such data set is discussed in the work of Koch. The work of Lindblad 
contains a related discussion for the fluid case. For the case of a Newtonian elastic 
body, we construct, for small values of Newton's constant, initial data satisfying 
suitable compatibility conditions near data for a relaxed elastic body in the ab- 
sence of gravity. It should be noted that due to the presence of the Newtonian 
gravitational potcintial this problem is non-local. 

The situation considered in this paper has several interesting generalizations. If 
we consider a large static, self-gravitating body, a relaxed state may not exist. Re- 
stricting to the spherically symmetric case, solutions to the field equations can be 
found by ODE techniques [28]. If we perturb the data slightly, there should exist so- 
lutions with small oscillations around the equilibrium configuration. Our approach 
allows one to consider also such more complicated problems; the essential difficulty 
is the construction of the initial data satisfying the compatibility conditions and 
not the time evolution. However, the results in this paper do no immediately cover 
these more general situations. In particular, a proof of a suitable version of the 
well-posedness result for the Cauchy problem, in these situations, requires a more 
general treatment of potential theory in Riemannian manifolds. This problem will 
be considered in a separate paper. 

Given a solution u to the Cauchy problem for the Newtonian elastic body which 
exists for for times t S [0, T], one may conclude from the main result of this paper, 
see theorem 4.3, that by choosing initial data sufficiently close to that of u, the time 
of existence of the resulting solution can be arbitrarily close to T. In particular, 
given a a static solution, there exist nearby data such that the corresponding solu- 
tion to the Cauchy problem has a time of existence not less than any given time T. 
The proof of existence of static self-gravitating bodies in Newtonian theory [6] could 
be used together with a generalization of theorem 4.3 that allows for more general 
initial data, as alluded to above, as well as a generalization of the construction of 
initial data to show that shows that nearby data define solutions which exist up to 
some time T. 

The treatment of two fluid balls in steady rotation due to Lichtcnstcin, which 
was mentioned above, has been generalized to the case of elastic bodies, see [8]. 
Using these solutions in the manner just described, it is possible to obtain classes 
of solutions of the two body problem in any prescribed flnite time interval. We 
leave these problems for later investigations. 
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Global existence for small data is an important question. For unbounded fluids 
in 3-dimensions, it is now known that shocks form for a large class of initial data 
[31, 11], and it is widely believed that for most equations of state shocks will always 
form for arbitrary small perturbation of constant density state. On the other hand, 
there do exists small data global existence results for unbounded clastic bodies 
provided the material satisfies certain additional conditions, cf. [32, 14]. For elastic 
bodies of finite extent, it appears to be an open question whether there exists 
solutions that are global in time or under what conditions shocks form. 

Overview of this paper. The paper is organized as follows. Section 2 describes 
Newtonian elasticity, sets up the basic equations and gives the conditions we impose 
on the material. We derive the equations in the material frame and in spacetime 
from a variational principle. Section 3 deals with the problem of finding solutions 
to the compatibility conditions needed for the proof of local existence. Due to the 
non-local terms in the equations, these condition imply conditions on the Cauchy 
data on the whole initial surface. To find initial data satisfying these conditions, we 
make use of some results from potential theory. These arc; developed in section 3.1. 
Further, the Poisson equation must be studied in the material frame, see section 
3.2. The results concerning the linearized elasticity operator which are needed can 
be found in section 3.3. Section 4 generalizes Koclrs theorem and proves our main 
theorem. Appendix A contains some background material for the function spaces 
used in this paper. 



2. Newtonian elasticity 

2.1. Kinematics. The body B is an open, connected, and bounded set with a C°° 
boundary in Euclidean space Mg. We refer to Mg as the extended body. We consider 
configurations, i.e. maps f -.M.^ ^ B, and deformations ^ : S — >■ with 

/o<^ = ide. (2.1) 

Thus, the physical body is the domain in space M| given by f^^{B) = <t){B). 

We shall make use of the extension ^ of to a map Rg — > K| and let {X'^)A=i,2,i 
and (2:')i=i.2,3 denote global Cartesian coordinates on the body Mg and the config- 
uration M5 spaces, respectively. Since we shall consider the Newtonian dynamics 
of a body, we let /, cf) depend on time, denoted by t. Equation (2.1) gives 

/^(i,<^(f,X)) = X^ inB and cj^^t, f{t,x)) = in f-\B). (2.2) 

Writing a;** = (f, x*), we introduce f"^^ = d^f"^ and (jj'^A = dA4''' ■ In particular, 
f\ = dtf^. We have 

/a/^ = (5^ and /^fc/A = ^''A (2.3) 
where these expressions are defined. This implies 

B^P'-A, and |^ = -/^J^. (2.4) 



Let 



, 1 in/-i(S) 
inK|\/-i(S) ^^■^> 
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be the indicator function of the support of the physical body. Using the above 
identities, one may calculate the variation of X/-i(S) '^ith respect to f^, 

^-^^ = <t>' Ad.Xf-HB) (2.6) 

Let H^^ = f'^,if^,j6'^ and define Hab by HabH^'^ = ^a^- Then Hab = 
4>^A<l>'BSij. We have 

J) IjAB 

1^ = 2S''-S(^c.f''K. (2.7) 

Differentiating (2.3) gives the usual formulas for the derivative of the inverse, 

d^,(l,\ = -<f>'Bd^f''k(l)''A , (2.8) 

dsf^" = -f^,dBclPcfi ■ (2.9) 

We let v^dfj, = dt + v^di be defined by v^f^^ = 0. This determines the vector 
field v^{x) on uniquely in terms of /. The velocity field describes the 

trajectories of material particles. From the relation v^dfj,f^ = 0, we get 

On the other hand, time differentiating (2.2) gives 

v\t,x) = {M'){tJ{t,x)). 

Thus, we have 

d^<l,''it,x) = dtiv\t,<j,it,x))) 

= {dtv''){t,cl){t,X))+dev''{t,<l){t,X))dt<l)^{t,X) 
= {v^^d^v''){t,<P{t,X)). 



Further, 

which shows that 



v'^v'^d^d^f^ = -v^^d^v^f^k, (2.10) 



HABV^^v'^d^d^f = -V^'d^v'^SmnrA- (2.11) 

The body B carries a reference volume element defined by a 3-form Vabc on B. 
The number density n is defined, cf. [3, eq. (3.2)], by 

f^rfofkVABcU[x)) = n{x)eiju{x) 

where is the volume clement of the Euclidean metric 5ij on M^. Since we are 
considering the Newtonian case with Euclidean geometry, we have simply 

n = det Df. 

The mass density is 

p = nm. (2.12) 

where m, is the specific mass of the particles, i.e. mV is the mass density^ for the 
material in its natural state, see [3] . Sec also [7] where the specific mass is denoted 
Pq. We have the relation 

QjA^=n<l>'A. (2.13) 



^Our techniques allow for m to be a function on the body, m = m{X), but for simplicity, we 
will assume that m is constant. 
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The following equations for n follow from the above definitions, using (2.13), 

d^inv'')=0, (2.14) 

and 

dk{n<f>''A) = 0. (2.15) 
The elastic material is described by the stored energy function 

The various forms of the stress tensor^ are defined from the stored energy function 



via 



We have 



de_ 



TAB = +'^-Q^AB ' '^ii = nf^if^jTAB 



de _i 



A. 



n TA 



de A 

- -Ti. 



d{(j)'A) 

The elasticity tensor Li^k^ is defined by 

It follows from this definition and the assumptions above that the elasticity tensor 
has the symmetries 

jJAjB _ j^jBiA _ j^iABj _ j^AijB ^2 17^ 

where 

Clearly, we also have 

2.2. Variational formulation. We derive the field equations for a self-gravitating 
elastic body from the elastic action, supplemented by a term giving Newton's force 

law and Newton's law of gravitation. The Lagrange density is of the form Aeoi23 
where €0123 is the 4- volume element on spacetime M x in Cartesian coordinates, 
and 

^ j^grav _|_ j^pot _|_ j^kin _|_ j^elast 

where 

xgrav l^^l 

~ SttG ' 

A elast 



■^The sign of the stress tensor here is the opposite of that of (Tab defined in [3, eq. (3.5)], but 
agrees with the usage in [7, equation (4.2)]. 
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with Xf-^B) given by (2.5), 

|VC/|2 = d^UdjU5'\ 

and 

= v^v^Sij. 

2.2.1. Eulerian picture. The action in the Eulerian picture then takes the form 

with A ~ A{U,diU, f^, f^o, f^i)- The Euler-Lagrange equations are of the form 
£a = 0, £u = with 

We have 

p AC/ 

Next, we consider the Euler-Lagrange terms generated by variations with respect 
to the configuration f-^. A calculation using (2.15) shows that the factor pXf-i{B) 
gives no contribution to the Euler-Lagrange equations, and hence the kinetic term 
in the action gives 

Ckin _ ^ d^kin _ 

which after some calculations, using (2.14) and (2.15), gives 

-fl" = -pv''d^v"'6mnrAXf-^(B) 

use (2.11) 

= pHABV>'v''f\f-^B)- 

The elastic term gives, using (2.15) and (2.6), 

\ elast n A elast 
_celast _ Q ±_ 

^ 'd{.f^,) d{.f^) 

= -9»(ta*X/-i(B))- 

In view of [3, lemma 2.2], we have that the divergence di{TA^X}-^{B)) is integrable 
only if the zero traction boundary condition 

'^^■'"■i 1 9/ -1(B) = 

holds, in which case the identity 

di{TA'Xf--^(B)) = diTA'Xf-^{B) 

is valid. 
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Finally, the potential term gives, using (2.15) and (2.6), 

= p4>'AdiUXf-i{B)- 

Adding the terms, we have 

-£a = pHABV^v^d^d^f^ - diT/ + p<f>\diU in f-\B), (2.18) 
subject to the boundary condition 

TA^nj\ j-i(g) = 0. 

Hence, we find that the Euler-Lagrange equations £a = 0,£^u = ^ siie equivalent to 
the system 

-pHABV^v'^d^d^f'' + d,TA' = P^'A^^U, in f-\B) (2.19a) 

AU = ATTGpXf-i(B), (2.19b) 

7-A*«j|9^-i(g) = (2.19c) 

where rij is the unit outward pointing normal to df~^{B). These are the field equa- 
tions for a dynamical elastic body in Newtonian gravity displayed in the Eulerian 
frame. Wc note that, by using (2.11) and multiplying by f^i, equations (2.19a) 
and (2.19c) take the form 

pvi'd^Vi + djTi^ = pdiU in /-1(B) and T^nj\g^_,^^^ = 0. (2.20) 
2.2.2. Material frame. Similarly, the action in the material frame is given by 
£™«* = J (j)*{Aeoi23)dX°dX^dX^dX^ 

= j JN^''*dX^dVABcdX^dX^dX'^ 

where 

J = {(t,*n)-^ = det{dA(t>') 
is the Jacobian of (j), and A"*"* = A""''(f/. <p'\<p\^^,<if a) is given by the relation 

j^mat ^ jlZ^ + {rnU + \mv'^ - e)xB (2.21) 
SttG 2 

where xb is the indicator function of the support of the body and C7, e are defined 
along the lines of [3] . In particular, U = U o (j), and 

|Vf7|| = H^'^dAUdBU 
is the puUback to Mg of |Vf7|^ where we recall that 

H^B ^fA,fB,gij ^.A ^ 

Similarly, e = e o ^ so that e = e{f^i, H^^). We note also that 

det(if^^) = det(/^0' = J^- 



8 



L. ANDERSSON, T.A. OLIYNYK, AND B. SCHMIDT 



The Euler-Lagrange equations in material frame are £i = 0, Eg = where 

Q£rnat Q £rnat 



Q£rnat Q£rnat 



'd{dAU) d{u) • 

This gives the system of equations 

-md'^cj)' + dA{f'^) = m5'^ ffdAlJ in B, (2.22a) 
AhU = AirGJ-^rnxB in M|, (2.22b) 
i^Af'-^loB = (2.22c) 
where va is the unit outward pointing normal to dB, 

f^^=<5*'=J(/^,rfe^)o<^ 
is the Piola transform of r,-' , and 

H = HAsdX'^dX^ = (j}*{6ijdx'dx^) 

is the pull back of the Euclidean metric Sijdx^dx^ under the map (f>. Observe that 
in (2.22b), we could also have used the notation p = J^^m since = no (p. 

Rescaling time and the Newtonian potential, we can write the evolution equa- 
tions (2.22a)-(2.22c) in the form 

-d^(l>' + dAf^' - X^S'^ffdAU = inB, (2.23a) 

AhU = J-'^mxB inM|, (2.23b) 
i^Af^'\oB = (2.23c) 

where 

We note that the Poisson equation (2.23b) can be written out more explicitly as 

dAiJE^^'dBU) ^mxB- (2.24) 
It follows from the symmetry properties (2.17) of the elasticity tensor that 

U^j^ = (2.25) 

^ ddB(l>> 

satisfies 

liAjB ^ pOiB ^ liABj ^ lAijB ^2.26) 

where 

Since we will be working in the material representation for the remainder of the 
article, we will drop the B from the body space and denote it simply by M^. 
For latter use, we define the following nonlinear functionals: 

E\cl>) = dA{f'^{dcl>)), and Ei{cj>) = tvaByAf'^{dcl>) (2.27) 

where dcf) = (9^^*). 
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2.3. Constitutive conditions. We make the following assumptions on the elastic 
material: 

(1) f^^ is a smooth function of its arguments {d(p) in the neighborhood of the 
identity map 

(2) the identity map is an equilibrium solution to (2.23a)-(2.23c) for A = 0, 
i.e. 

and 

(3) 



„iA B _ fiA B 
il j — J-j j 



8<p=8if!o 



satisfies the following properties: 

(a) there exists an a; > such that 

a'^^''{X)^,UVjVB>u>\e\v\' 

for all ^,77 G R^, and 

(b) there exists a 7 > such that 

la'^^^'diAajB > 7kl' 

for all a = (cTja) G M3X3 with ats = c^Bi- 

We remark that condition (2) above is satisfied for a stored energy function which 
has a minimum at some reference configuration. See [3, equations (3.20)-(3.22)]. 

3. Construction of initial data 

In order to prove the existence of dynamical solutions to the evolution equations 
(2.23a)-(2.23c), we first need to construct initial data that satisfy the compatibility 
conditions to a sufficiently high order. The compatability conditions are defined as 
follows. 

Definition 3.1. Fixing s > 3/2 + 1, we say that the initial data 

{cl>'\t=o,dt\t=ocl>') - (cl>lcl>\) G W'+'^\B,W') X W'^^(B,R^) 
satisfies the compatibility conditions to order r (0 < r < s) if there exists maps 
(f)leW'+^-'^'^{B,M.^) e = 2,3,...,r 

that satisfy 

df-^ {-dW + dAf^' - X^S'^ffdAU) 1,^0 = in B, 



p3 



for £ = 0, 1, 2, ... r where, after formally differentiating, we set df\t=o(i 
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3.1. Potential theory. Before we can solve the problem of existence of initial data 
satisfying the compatibility conditions and the time evolution of this data, we first 
need to develop some potential theory. To begin, we set B+ = B, B- = M.^\B, and 
let 

E(X, Y) = ; 

denote the Newton potential so that AE = S. We also let S and V denote the 
single and double layer potentials 



S[f]{X)= [ E{X,Y)f{Y)da{Y) X ^ dB 



and 



V[f]{X)= ( ^E{X,Y)f(Y)da{Y) X OB, 
Job C;^y 

respectively, where da is the induced surface measure on dB, and v is the outward 
pointing normal. Restricting to dB, we have for X G dB that 

S[f]\^^{X) = J^^E{X,Y)f{Y)da{Y) and V[f]\^^^{X) = {±^1 + K)[f]{X) 
K[f]{X) = P.V. J^^ ^E{X, Y)f{Y) da{Y). 



where 



Further, 

where K* is the adjoint of K. We recall the following well known relations between 
the boundary value problems for A in B± and these potentials: 

(i) The solution to the Dirichlet problem 

Au = 0, tvgsu = V' 

on B± is given by 

u = V[f] 

where f solves 

(±i/ + if)[/] -V- 

(ii) The solution to the Neumann problem 

Au = 0, tVdB —u = ip 
ov 



on B± is given by 
where / solves 



u = S[f] + C 



2 

and C is an arbitrary constant. The solution exists if and only if J^g ip da 
0. 
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For the moment, we consider the interior Dirichlet and Neumann problems on 
B+ = B. The solution for the Dirichlet problem has the property that 

u e W''P{B) if troBU G B'-'^^P'P{dB), (3.1) 

while for the Neumann problem, we have 

ueW''P{B) if trae -^u e B'-'^ " ^P'P (dB). 

These results are classical, sec [2] . Defining the volume potential of a density / in 
Bto be 

V[f]{X) = A-\fxB){X) = f E{X,Y)f{Y)d'Y, 

Jb 

we abuse notation and say that V[/] G W'' p{B) if the restriction to B has this 
property, and similarly for the other potentials. Using the fact that dxAE{X, Y) = 
—dYAE{X,Y), we have 

dAV[f]{X)^- I dYAE{X,Y)f{Y)SY. 

JB 

This gives, after a partial integration, 

dAV[f] = V[dAf] - S[tvdBfA- (3-2) 
Let u be the solution of the Dirichlet problem 

Au = 0, tidBU = ip- 

By Green's theorem, wc have 

/ {Ayu{Y)E{X, Y)-u{Y)AyE{X, Y)) d^Y = 
Jb 



Since u solves the Dirichlet problem with boundary data ^, and 

AyE{X, Y) = S{X - Y), (3.3) 

this gives 

d 

V[u]\B=S[tioBg^u]+u. (3.4) 

Upon taking the limit from the interior at dB, we have 

(i/ + K)[u]\aB = S[tion£u] + u. 

Consider a metric qab on with covariant derivative V^. Let {ea}a=i,2 be a 
tangential frame on dB, and let hab denote the induced metric on dB with covariant 
derivative D. Let a vector field ^ be given. Decompose ^ into tangential and normal 
component at dB, 

Introduce a Gaussian foliation near dB. Then g takes the form 

gAEdX^dX"" = dr^ + habiv, r)dy-'dy'' 

where are coordinates on dB. Extending j/ in a neighborhood of dB using a 
Gaussian foliation, we have 

V^^^ = /i«*V„6 + (V0(z^,'^)- 
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The last term vanishes due to V^j^ = 0, which is vahd in a Gauss foliation. Then 
we have 

where Xab = (VoZ^, Cb) = \Cvhab is the second fundamental form. Let H = h°'''Xab 
denote the mean curvature of dB. Then we have 

Now specialize to the Euclidean case; let qab = 5ab be the Euclidean metric and 
let Oa for some fixed A. Then Va^"^ = 0, and we have 

DaP" = Hv^. (3.5) 

For X e S, we calculate using (3.3) and the divergence theorem that 

dAS[J]{X) = - I J{Y)P-DaE{X, Y) da{Y) 

- J^^ m {y, dy. )-^E{X, Y) da{Y) 

= / (D^P'^fiY) + P''DJ)E{X,Y)da{Y) 

JOB 

- [ f{Y){v,dy.)^E{x,y)dcj{Y). 

JdB 

Thus we have by the above result and (3.5) that 

dAS[f] = S[fHu^ + dlf]- V[fu^] . (3.6) 

Proposition 3.2. The operators S,^I + K have the mapping properties 

ll + K ■.B''-'^/P^P{da) S'^-i/P'f (5^), (3.7a) 

S (an) (an), (3.7b) 

for s = k — 1/p, k > 1, k an integer. The corresponding statements for the single 
and double layer potentials are 

V: S'=-VP.P(an) ^- W'='P(f2) (3.8a) 

S: S'^-^-i/f'f (an) ^"^='^(0) (3.8b) 

for k>l,k integer. 

Proof. We use induction to reduce the statement to the case A; = 1. This case 
follows from [15], see also [26]. Suppose then that we have proved the statement 
for A; — 1. To do the induction, assume 

/ g B'=-i/P'P(an). Let u be the solution 
to the Dirichlet problem with boundary data /. Then u G W'''P{fl), and setting 
= trgn^u, we have ip G B'^^^^^^p^p . Equations (3.4) and (3.6) give 

d,,v[f] = d^,s[^] + d,,u 

= S[iiHv' + af,V] - VliPv^] + ^^^u 

which by the induction assumption is in W^~^'P . It follows that V[f] G W^'P{^) 
and hence (i/ + K)f G B''~^/P'P{d^). This proves the statement for for T) and 
{^I + K) at regularity k. 
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Next, for S, we use the equation (3.6) for / G B'^~^~^/P'P{dQ). Using the 
induction assumption and the statement for D and (5/ + K) just proved, we have 

which gives S[f] e W'''P{dn), and S[f] G B'=-Vp.P(af2). This completes the induc- 
tion and the result follows. □ 

Example 3.3. Let / e W^'P{n). Then tionf e B^-^/p^p and V[d^..f] £ W^'P{n). 
Further, iS[traj2/i^'] solves a Dirichlet problem with boundary data 5[trao/i'*] G 
B^-^/P'P{dn), and hence <S[traQ/i^'] € W^^p. It follows, in view of (3.2), that 
(V/)xn G M^3,p(f^)_ 

Example 3.4. Let / e PF2'P(S). We have 

Since 5^./ G V^i'P(i3), we have from Example 3.3 that V[d^,f] G W^'P{B). Further, 
traa/ G B2-i/p,p(c)0) and hence 5[trao/i^*] e H^^.p^^^^ Therefore V[/] G W^'P{B). 

Proposition 3.5. Let k>l, and assume f G M^'='P(S). Then V[f] G W'=+2'P(B). 

Proof. The proof proceeds by induction, with base case k = 1. For this case, 
the statement follows by the argument in Example 3.3. Suppose wc have proved 
the statement for fc — 1. We will make use of the identity (3.2). By induction, 
V[d^,f] G W''+^^P{n). Further, tran/ G B'^-'^/p and hence by (3.7), 5'[trao/i/*] G 
Bk+i-i/p(^On). It follows by (3.1) that S[tianfv'] G W''+'^'P{B). This shows that 
dx'V[f] G W'+^'PiB) and hence V[/] G W'+^'PiB). □ 

Similar arguments combined with the mapping properties (A.7)-(A.8) of the 
Laplacian on the weighted Sobolev spaces can be used to establish the following 
proposition for the volume potential on jB_. 

Proposition 3.6. Let k > 1, -1 < 5 < 0, and assume f G Wgf^i^-)- Then 
V[f] G +''^(B_). 

3.2. The Poisson equation in the material frame. The next step in solving 
the problems of the existence of initial data satisfying the compatibility conditions 
and the time evolution of this data is to establish a number of smoothness properties 
for solutions to the Poisson equations in the material frame. We begin by defining 

the spaces 

yy;M,P(^3) ^ | ^ g w^^^ {M.\ V) \ u\b G I^'^+^^f (B) and u\b_ e W^+'^P{B-) } 

for 1 < p < 00, A: G Z and 5 € R. It is not difficult to verify that these spaces are 
complete with respect to the norm 

ll^ll>V^*'"'''(K3) = I|w|b||w^'=+-.p(B) + I|w|B-IIh/^*^+ = 'P(B_) + ll'"llw^'='f(R3)> (3-9) 

and hence Banach spaces. 

Theorem 3.7. Suppose 1 < p < 00, s G Z>o, s + 1 > 3/p, and — 1 < ^ < 0. Then 
there exist an open neighborhood O^+^'P c W^~^'^'P{B,R^) of tpo, and an analytic 
map 

such that 
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(i) U{(p) satisfies the Poisson equation (2.23b) on R^, 

(ii) for each (f) G O^+^.p^ the map cj) = ^o+IEb(^— '^olS) is a diffeomorphism 
onIR3 that satisfies 4>-il>o € w!|o'''(I^^. 1^^) C Ciio(IR^. 1^^) and4>-^-^o € 

(iii) [/ = [/((/>) o G W^^'^(IR.^) satisfies the Poisson equation AU = mpx4,{B) 
on R? and U{x) = o(|a;|'') as \x\ — >■ oo, and 

(iv) for fc e Z>i and k < s + 1, the derivative of U can he extended to act on 

W^''P(B), and moreover, the map 

0«+2,p 9 ^ ^ DU{(j)) G L(W^*='f (S)) 
is well defined and analytic?. 

Proof (i) Fix 1< p < 00, s + 1 > 3/p, and -1< 5 < 0. Given ^ G W^+^'^iB), we 
define 

= Vo+Eb(V'), <l>\ = dA^\ and (//i) = (</>^)-i. 
Since matrix inversion and the determinant both define analytic maps in a neigh- 
borhood of the identity, it follows from (A. 6), proposition 3.6 of [16], the continuity 
of extension and differentiation, the analyticity of continuous linear maps, and the 
property that the composition of analytic maps are again analytic that there exists 
a > such that the maps 

BR{W^+''P{B,m?)) 9 V ^ det((^^) - 1 G W^Xl'^i^^) (3.10) 

and 

Br{W^+^'%B,R?)) ^^^{ff- 6f) G W^1+^^^(R3,M3x3) (3.11) 
are well defined and analytic. Recalling that H^^ = f^S^^f^, we see from the 
same arguments that the map 

Br{W^+''P{B, R3)) 9 V ^ JH"^^ - e W^l|l'■^'(M^ Msxs) (3.12) 

is analytic where J = det{^\). Using the multiplication inequalities (A. 2) and 
(A. 6), we find that 

\\dA{JH^''dBU)\Jw'.,(s) < \mw^+..iB,M^)\\U\B\\w^+..,^s), (3.13) 
||aA(Ji?-'^5BC/)|e_||w;-(B_) < \mw^^'-W^)\\U\BJw!+^'^B.) (3-14) 

and 

\\dA{jH'^^dBU)\\^o^P^(^s) < ||V'|k2+.,P(B,R3)||[/||vi/|,P(K3). (3.15) 

From the analyticity of the map (3.12), the bilinear estimates (3.13)-(3.15), the 
analyticity of continuous bilinear maps, and the property that the composition of 
analytic maps are again analytic, it follows that the 

Br{W^+''P{B,M.^)) X W|'*'^'(R3) 9 (^,[7) ^ dAiJE^^'dBU) G W°!_Y(K') 

(3.16) 

is well defined and analytic. Together, (A. 7), and Propositions 3.5 and 3.6 imply 
that 

A-i(xs)eW|'^'^(R^), (3.17) 

and the Laplacian 

A : Ws'^'^iM.^) — s- W°:_Y(R^) (3.18) 
•^For a BaiiELch spELce X, L(X) denotes the set of continuous linear operators on X 
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is an isomorphism with inverse given by (A. 8). 
Prom (3.16), (3.17), and (3.18), we see that 

(V, C7) ^ (Oa (JH^^'dBU) - rnxs) 
is well defined and analytic. Evaluating F at V = gives 

F{0,U) = A-^(^AU-mxB), 

which shows that 

Do = mA-\xB) e Wl''P{R^) (3.19) 

satisfies 

F(0,C7o) = 0. (3.20) 
Also, b}' the linearity of F in its second argument and the invertibility of the 
Laplacian, it is clear that 

D2F{0, U)-SU = (A(5C7) = 6U. (3.21) 

Results (3.20) and (3.21) allow us to apply an analytic version of the implicit 
function theorem (see [13], theorem 15.3) to conclude the existence of a unique 
analytic map, shrinking R if necessary, 

tJ : Br{W'^+''P{B,M.^)) ^Wl''P{R^) (3.22) 

that satisfies 

U{0) = Uo, 

and 

F(V', U{i;)) =0 V V G Bii{W'^+''P{B, R^)). 
From the definition of F and the invertibility of A~^, it then follows that (/{tp) 
satisfies 

dAiJE^^dBUi^)) =mxB- (3.23) 

(ii) & (iii) Following Cantor [9] , we consider the following group of diffeomorphisms 
on R^ 

X>^'«(M3) :={(j):R^ ^R^\(j)-2pQ€ W^*'9(M^M^), and (p-^ - Vo G W^''«(K^R3) } 

where s > 3/g + 1 and ?? < 0. Fixing tp G Br{W'^+^'P{B,R^)), we get from (3.10) 
and (3.11) that 

Defining, 

U ■.= U{ip)o^-\ (3.24) 
we can apply Corollary 1.6 of [9] to gef* 

U e Ws'^iR"^). (3.25) 

A straightforward calculation using the chain rule and (3.23), (3.24), and (3.25) 
then shows that 

AC/ = mdet(i?(r'))x^(B), 

"^In [9], Cantor required that <5 < —3/2 because that was what he needed to prove the weighted 
multipUcation inequahty (A. 6). It is clear that his proofs are valid whenever the multiphcation 
inequaUty holds and Wg'^ C C^. Consequently, the only restriction on S is that <5 < 0. 
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while the faU off condition U{x) = o{\x\^) as \x\ — ?■ oo follows from the weighted 
Sobolev inequality (A. 5). 

(iv) To begin, we assume that k = 1 and observe that for 0^ e W^'P{B) and 



< ||^||w'i.P(B)||t^|B|lw^5"+2.f(B) + II^IIi^1'''(B)||£''|B-|Im^^-+2.p(b_) 

< ||^||wi.P(B)l|t^llw2.^.''(R3) (3.26) 

where in deriving the result wc have used property (A. 9) of the extension operator 
Eg, the multiplication inequalities (A. 2) and (A. 6), and the assumption 1 + s > 3/p. 
Letting 

ttAbc _ dJH^^ 

we get, using the estimate (3.26) and the same arguments as above, that for R 
small enough the map 

Ga : BR{W^+''P{B,m^)) X W^'P{B,W^) x W|'"'^'(R3) x W^'P{m^) W^^f{W^) 

: e\ u, V) ^ jH^'^dBV + Hf''^EB{dce^)dB{u) 

is analytic. From the continuity of differentiation and the trace map, we then have 
that the map 

G : Br{W^+''P{B,R^)) X W^'P{B,M.^) x wf^'^iW^) x W^^'^(R3) W^^fiW^) 
defined by 

G{r,9\U,V) = dAG^{r,0'B,U,V) 

is analytic. Taking Uq as defined by (3.19), a straightforward calculation and the 
invertibility of the Laplacian show that 

G(0,0,C7o,0) = and DiG{0,0,Uo,0) ■ SV = SV. 

Therefore, wc can again apply the analytic version of the implicit function theo- 
rem (see [13], theorem 15.3) to conclude the existence of a unique analytic map, 
shrinking R if necessary, 

V : Br{W'^+''P{B,R^)) X W^'P{B,R^) x (C/q + ^^(W^'''^)) W^'^iR^) (3.27) 

that satisfies 

y(o,o,[/o,o) = o, 

and 

G{r,0's,U,V{r,0's,U))=0 

for all {il}\e\U) e Bii{W^+'-P{B,R^)) x W^-p{B,R^) x {Uo + Br{wI"'P{R^))). 

From the construction of G, and the uniqueness of the maps (3.22) and (3.27), 
it is not difficult to verify that 

Vi^\dA6^p\Uii:)) = DUi^P)-S^P V(V',(5V') e BR{W^+''''{B,M.^))xW^+''PiB,R^). 

From this and the density of W^+^'^'(i3) in W^'P{B), it follows that the derivative 
of U can be extended act on W^'P{B), and moreover, that the map 

Br{W'^+''P{B, R^)) 9 V i-^- DU{tp) e L{W^'P{B)) (3.28) 
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is well defined and analytic. By (3.22) above, wc also have that the map 

BniW^+'^'PiB, R3)) 9 V ^ DUiij) e LiW^+^'PiB)) (3.29) 

is well defined and analytic. Together, the maps (3.28)-(3.29) and interpolation 
imply that the map 

is well defined and analytic for fc G Z and l<fc<s + 2. □ 
CoroUciry 3.8. The map 

A : o«+2,P — ^ VF^+i'P(B,M^) 

defined by 

is analytic. Moreover, for k G Z>i and k < s + 1, the derivative of A can be 
extended act on W'''P{B), and the map 

O'+^-P DA{^) G L{W'''P{B),W''-^'P{B)) 

is well defined and analytic. 

Proof. This follows directly from theorem 3.7, the multiplication inequality (A. 6), 
proposition 3.6 of [16], and the fact that compositions of analytic maps are again 
analytic. □ 

3.3. The linearized elasticity operator. We define the operator linearized elas- 
ticity and boundary operators by 

A(</))* :=aB(a*^,-^9c<^) (3.30) 

and 

Aa{<j>y := :/Ba*^/5i5<^, (3.31) 

respectively. We also define 

Y''P = W''P{B,R^) X B'+'^-'^/P'P{dB,R^), 
and for each G W'+'^'P{B, R^), 

y;'P = { (h, t) G r^'f I Ci(5, t) = 0, C2((/., h,t) = o} 

where 



Ci(6,t)= / 5+ / t, (3.32) 

JB JdB 

and 

C2{(l),b,t)= bx(l)+ tx(j). (3.33) 

JB JdB 

Here, we are using the notation 

(5 X (l>y = (3.34) 



f b = [ bd^X, (3.35) 
Jb Jb 
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and 

/ t = I tda. (3.36) 

JdB JdB 

For later use, we recall the following theorem from [24] concerning the surjectivity 

of the linearized elasticity operator. 

Theorem 3.9. [[24], theorem 1.11, Section 6.1] Suppose 1 < p < oo, s > Q, and 
P : — > Y^'^ is any projection map. Then the map 

A : W^+^'P{B,R^)^Y;f : ^ P(A(0), 

is surjective and 

ker A = {a + bxX\a = {a'), b = {V) e M^} 
where {b x Xy = e'jAVX"^. 
Remark 3.10. Letting 

denote the center of B, a short calculation using the change of coordinates X^ = 
Xi -X^ anAB = B-X shows that 

/ X^d^X= [ X^-X^d^X. 
Jb Jb 

Therefore, we can always arrange that 

X^ d^X = (3.37) 



L 



IB 

by translating the domain B. For the remainder of this article, we will always 
assume that the condition (3.37) holds. 

Next, for p > 3/2, we define the spaces 

x^'P = ^(l)eW''P{B,M.^) y"(^ = o}, 

and observe, using Sobolev's and Holder's inequalities, that 

(I) X tjj 

Jb 

^ " ' W2'f(B,K3)||V||LP(B,I 
< [[0||vy2+».P(B,M3) ||'^i'||iy»>f (B,IR3) 

from which the continuity of the bilinear map 

B ■ X X* — )■ : ((^, V-) ^ / X -0 



< ll"?^ X V'l|z,i(B,M3) < ||?^||loo(B,r3)[[V'||li(B,K3) 

< 



follows. Setting 

we define for p > 3/2 the following spaces 

W'P = {il)G X''P [ : ker A n X'+'^'P is an isomorphism }. 

Lemma 3.11. Vo € W'P for all s>0 and 3/2 < p < oo. 
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Proof. By (3.37), we have that 

/ tpid^X= [ X'd^X = 0, 

J Be Jb 

and 



/ a + bxX= [ a + bx [ X = Vol(e) 
Jb Jb Jb 

for all a, 6 G K^. Consequently, 

V'o e X'^P, (3.38) 

and 

kerAnX''+^'P = {6xX|6eM^} (3.39) 

by theorem 3.9. Next, 

B^^{b xX)= [ {bxX)xtpo= [ {bxX)xX 
Jb Jb 

= / iX-b)X-b\Xf, 
Jb 

which, after taking the innerproduct with a e M^, yields 

b-B^^{bxX)= [ {X- b){X -a) -a- b\Xf. (3.40) 
Jb 

The Cauchy-Schwartz inequality shows that 

{X -bf -\bf\Xf <0 (3.41) 

and 

{X -bf -\b\^\X\^ = VXGB^b = 0. (3.42) 
Combining (3.40)-(3.42), we arrive at 

b-B^„{bxX) = 0^b = 0. (3.43) 

By way of contradiction, suppose that the map 

ker A n — )■ (3.44) 

is not surjective. Then the image B^^ (ker A n X*+^'P) is contained in a two dimen- 
sional subspace, and therefore, there exists a non zero a G such that 

a-B^„{bxX)=0 V6gK3. 

But this is impossible by (3.43), and hence the map (3.44) is surjective. Since 
dim ker A n X^^'^'P = 3, the map (3.44) must, in fact, be an isomorphism. □ 

3.4. Existence of initial data satisfying the compatibility conditions. 

Lemma 3.12. Suppose that J^o; 3^i) ■ • ■ ) 3^r) (ind Z are Banach spaces with contin- 
uous (linear) embeddings 

H,j -yi — >yj i,j G{0,l,...,r}, i< j, 

U is open, and F G {U, Z) . Then the map defined by 



Fr{y^i, yi,...,yr) ■■= ^0^^^^^^^ '^^^^^ ^^^^ " ^jAVi) 



dt 

isinC^{i^^,{U)xYY'.^^y^,Z). 
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Proof. Since U £ Yr is open, it follow from the continuity of the map Lo,r that 
toliU) C lo is open. Next, fix j/o € '-o r(^) Vj ^ 3 — 1, 2, . . . , r. Then 

the continuity of the maps Lj^r '-Yj and (.o,r(yo) € U guarantees the existence 

of a 5 > such that 

r 

c{t) =Y,t' hAVj) &U ^t€{-S,5). 

3=0 

Clearly, this implies that c € C°°{{—5,6),U), and hence, that the map 

C(^) X n 9 {yo, .■.,yr)^ ^1=0^^""^*^^ ^ ^ 

is well defined and continuously differentiable. □ 

Proposition 3.13. Suppose 3 < p < oo, s € Z>o, and O^+'^'P C W^+^'P{B,R^) is 
the open neighborhood o/V'o from theorem 3.7. Then the maps (see (2.27)^ 

Eg : 0«+2.P — ). B*+i-Vp.p(5S) 

are C°°. 

Proof. First, we recall that, by assumption f^^, is a smooth function of its argu- 
ments in the neighborhood of the identity map tpQ. Since p > 3 and s > 0, 
we have that s + 1 > 3/p, and it follows from theorem 1, Section 5.5.2, of [30], and 
the continuity of differentiation (cf. (A.l)) that the map C)-'+2.p 9 (+, 1 — y r*'^ G 
W^~^^'P{B,R^) is C°°. The proof then follows directly from the continuity of the 
trace map (A. 3). □ 

Defining 

QS+2,p _ Qs+2,p Pi j^s+2,p 

we have that 

V-o e 0"+2,p 

by (3.38). We also define^ 
F : O'+^'P X X''P X K ^ y^'f X . (0o,</.2, A) ^ (7-(0o,<^2, A),B((Ao,V'o)) 
where 

/•((^o, '^'2, A) = (i;(0o) + a2A(0o) - (/)2, i;a(<Ao)) , 

and we let 

p : O'+'^'P L{Y'-P, Y'^P) (3.45) 
denote any C°° map for which P^^ coincides with the projection operator from 
theorem 3.9. Furthermore, we assume that for each ^0 € O^+^'P, the linear operator 

is an isomorphism and 

y«'P = y|'P©kerP((/.o). (3.47) 
The existence of a map (3.45) satisfying (3.46) and (3.47) can be found in [20]. 



^For two Banach spaces X and Y , L(X, Y) denotes the set of continous linear maps from X 
to Y. 
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For r e Z>o, we define 

j=o •^^ 

and 

-Po(<?^o, <^2, A) = F((?io, <?i2, A), 

.^3, A) = ^ (l>s, A), .^i,2(t), A), AS(Vo, (t>i)) , 

If cP~ 

-PV+2(</>r,0r+2, A) = - ( j;+2(0r+2)'?^'r+4, A), — ^_^B(0o (0 > 00,2 (*) ) 
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t=0 



(r > 0) 



where 



j;(0^,(/)r+2,A) = — ^ ^J^{(j)r{t),(j)r,2{i)^^)- 



t=0 



Remark 3.14. Under the identification 

<^r = (ar</>)|,^o, 

it follows from the definition of 4'r{t) and (j)r,2{t) above that 



t=o 



yv^ii. „ and -— 



t=0 



^2(i) = A(a[+V)L=o- 



Moreover, under this identification, we have that 

7-(<^o,</'2,A) =0 

if and only if 

{dAf'^{d4>) - \H'^ffdAU{<j>) - dU') l*=o = and {vAf'^{d<t>)\oB) |t=o = 0. 

Also, by repeatedly differentiating the equations of motion (2.23a)-(2.23c), it is not 
difficult to see that 



(3.48) 



if and only if 



dl{dAT'^idcb)-X^P'ffdAUi<P)~d^cj)%^o = and a[ (z^^f^^ (90)1 9^3) |t=o = 0. 

This shows that solving (3.48) for r = 0,1,...,^ will produce initial data that 
satisfies the compatibility conditions to order i. 

In order to use the implicit function theorem to solve the equations (3.48), we 
need to introduce the following maps which are a closely related to the J> and Fj- 
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maps introduced above: 
^(</.o,</'2,A) = P(,^oWo,02,A), 

Gi(</.i,03,A) = i(gi(<^i,03,A),AS(^o,0i)) 



Gr+2(^r+2> 0r+4, A) 
1 /d^+^ 



A V 



^gr+2{4>r+2{t), 0r+4, A), ^ |^^^B(</)o(t), 00,2^)) , (r > 0) 



and 



Gr(0r>0r+l,</'r+2,A) 



where 



/ Go(0o><^2,A) \ 
Gi(</.i,03,A) 

\ Gr(0r,^r-+2, A) / 
5r(0r)'/>r+2,A) = — 0r,2(t), A). 



Proposition 3.15. Suppose s G Z>r and 3 < p < oo. T/ien i/ie maps 

r+2 



and 



are C^. 



r+2 



Proof. First we note that the maps Fq and Go are C°° which follows from Corollary 
3.8, proposition 3.13, and the smoothness of the map (3.45). The proof then follows 
immediately from lemma 3.12 and the definition of and Gr- O 



e{(l,,\) = {E{4>) + \^K{<j>),Ee{<i>)), 



Introducing 
we get that 

i^o(<^o, 02, A) = {S{4>o, A) - ((/.2, 0), B(V'o, </>o)) , 
and it follows easily from the definition of the Fr maps above that 

Fi{cl>^,cl>3, A) = {D^£i(l>o, A) • - ih, 0), B(Vo, , 

and 

^^2(02, cl>4, A) = {D^£{cl>o, A) • 02 - (04, 0) + XDl£{<f>o, A) • (0i, 0i), B(0o, 02)) • 
Proceeding inductively, we obtain for r > 1, 

JV(0^,0^+2,A) = (£>^£(0o,A) -0^ - (0^+2,O),B(0o,0r)) +AF^(0^+i,A) (3.49) 
where the map Fr is . Setting, 



Ki<Pr) = ^ 



t=0 
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the product rule shows that 

gr{<Pr,<t>r+2,X) = ^ ('']Fr-j{<Pr-j)M4>jAj+2,X). (3.50) 
3=0 

This and (3.49), in turn, give 

Gr{<Pr, <Pr+2, A) = (P(.^o) [Dc^S{(I>0, A)-</'r-(</'r+2, 0)] , B{(t>o, .^^))+AGr(0^+i, A) r > 1 

(3.51) 

where Gr is C^. 
Letting 

= (^o,0,...,0), 

it follows directly from (3.51) that the derivative of G^. evaluated at (0^, (j)r+i,(j)r+2, A) = 
(V^, 0,0,0) is 

D^^Gr{tpv,0,0,0)-S(l)r=A-5(l)^ (3.52) 

where 



fCi^o) MiiJo) 

C{i;o) M(^o) 

£(^-0) 

C{iPo) 



\ 





V •■• C{tPo)/ 

C{i;o)-6i; = (P(^^o)A(5^/;,B(V'o,^^)), 



(3.53) 



(3.54) 



and 



A^(7Ao)-'5V' = (-F(V'o)('5V',0),0). 

We note that in deriving (3.52)- (3.53), we have used 

£{tpo,0) = and D^£{^l)o, 0) ■ 5^p = 

Although our main objective is to solve the equations (3.48), we first solve Gr = 
and later show that this implies that compatibility conditions arc satisfied to order 
r. To solve G^. = 0, we use the implicit function theorem. The proof we present is 
modeled on the existence proof for static, self-gravitating elastic bodies presented 
in [7]. 

Proposition 3.16. There exists an Aq > 0, open neighborhoods Nr+i C X^'^^~'^'P 

and Mr+2 C X'^^'^'P both containing 0, and maps 

: 7V;+ixA^r+2x(-Ao,Ao) — ^ x^+^-^'-p : (0^+1, ^r+2, A) i — s- $j((/)r+i, (/)r+2, A), 
for J = 0, 1, . . . , r, such that 

$o(0,0,0) = V^o, $j(0,0,0) = 0, j = l,2,...,r 

and 

Gr{^r{<t>r+l,(}>r+2),(t>r+l,<f>r+2,>) =0 V((?!>r+1 , (/'r+2, A) G A/'r+l xA/'r+2 X (-Aq, Aq) 

where *r = (^0,^1, . . • ,^r). 

Proof. We begin by verifying the the linear operator (3.53) is an isomorphism. 
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Lemma 3.17. For s G Z>r and 3 < p < oo, the map 

r r 
j=0 j=0 

is a linear isomorphism. 

Proof. By lemma 3.11, tpo G W'^ and so it follows from the definition of jC{tpo) (see 
(3.54)) and theorem 3.9 that C{ipo) ■ X'+'^'^'P — > y*-J> xM^ jg an isomorphism for 
j = 0,1, . . .r. The proof now follows immediately from upper triangular structure 
of the map A (see (3.53)). □ 

Recalling that £{tl>o) = 0, it is clear from the definition of T that 

jr(V;o,0,0) = 0, 
and hence, by the antisymmetry of the map B, that 

Go(^o,0,0) = 0. 
Furthermore, it clear from (3.51) that 

G,(V',-,0,0)=0 j = l,2,...,r, 
which, in turn, shows that 

Gr(V'r,0,0,0) =0. (3.55) 
The proof of the proposition now follows proposition 3.15, lemma 3.17, and the 
implicit function theorem. □ 

Theorem 3.18. The m,aps $j j = 0, 1, . . . , r from, proposition 3.16 satisfy 

J'j(*j(4-+l,4-+2),$j+2(0r+l,0r+2),A) =0 j = 0, 1, . . . , r ~ 2, 
-^r-l(*r-l(0r+l, 0r+2), (Ar+1, A) = aud 7"^ (*r(</>r+l, ?!'r+2), ?f'r+2, A) = 
for all {(f)r+l,cj)r+2, A) G Nr+l X Mr+2 X ("Aq, Aq). 

Proof. It is shown in [7] that 

Ci{£{(l),X)) =0 and C2 (0, f A)) = (3.56) 
are automatically satisfied for all cf) G O-'+^'P and — Aq < A < Aq. Setting 

t=o 



£j{<j>j,X) = ^ f{<l>r{t), A) j = 0, 1, . . . , r, 



the formulas 



Ci(f,(0^.,A)) =0 and ^ [^y2(</'fc,^i-fe(0,_fc, A)) = j = 0,l,...,r 

follow from differentiating (3.56). In particular, this implies that the maps $j from 
theorem 3.16 satisfy 

Ci(£,(*,(.^,+i,</.,+2,A)),A)) =0 for j = 0,l,...,r (3.57) 
and for J = 0, 1, . . . , r — 2, 

51 (i )'^2(*fe(^r+l,'/>r+2,A),fj_fe(*j_A;((/)r+l,(/)r+2,A),A)) =0, (3.58) 



fe=0 



for all ((^r+l, <^r+2, A) e Nr+\ X Nr+1 X (-Aq, Aq) 
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Next, we note that 

Ci((<^,, 0)) =0, y(t>r& (3.59) 

and 

dP 



^^^i3(0,(i),02,,(t)) = 0^ ^[^^C2(</'iW,(</>2,rW,0)) = 



^ E (i)c2{h,{^j+2-k,0)) = 0. (3.60) 
fe=o ^'^^ 

If we define 

$r+l('/'r+l, </'r+2, A) = and ^r+2{4>r+l, (pr+2, X) = (pr+2, 

then it follows from (3.59), (3.60), and the definition of Gr that the maps satisfy 
Ci($j(0,+i,0,+2,A),O) =0 j = 0,l,...,r + 2, (3.61) 

and 

(''j^C2{^ki(f>r+l,<f>r+2, X), {^j+2-k{<t>r+l,(l>r+2, =0 i = 0,l,...,r 

k=0 ^ ^ 

(3.62) 

for all (0r+l, (/>r+2, A) G A/'r+l X 7Vr+2 X (-Ao, Ao). 

Fixing ((^,.+1, (/)r+2, A) G Afr+i x AAr+2 X (-Ao,Ao), the identities (3.57), (3.58), 
(3.61), (3.62) and the definition of the maps J-j show that 

Ci(J,(*,(<^,+i,<^,+2,A),$,+2((/..+i,.^.+2,A),A)) =0, (3.63) 

and 

V] ( 'l]C2{^k{(f>r+l,(f>r+2, X) , Tj-k{^ j-k{4>r+\, 4>r+2, A), <^ j+2-k{4>r+\, 4>r+2, A), A)) = 
fe=0 ^'^^ 

(3.64) 

for < j < r, while 

^o(^o(0r+l,<?^r+2,A),$2(<?^r+l,'?^r+2,A), A) = 

follows from theorem 3.16, or equivalently, by the definition of Qo, and 

P($0('?^fc+1, '?^r+2, A)) J"o(*o('?ir+l, 'Ar+2, A), $2(9^r+l, 0r+2, A), A) = 0. (3.65) 

But, we also have that 

J"o(*o(<?^r+l,<?ir+2,A),$2(<?^r+l,0r+2,A),A) G ^$;^0^_^j,^^_^,,A) 

by (3.63)-(3.64), and thus, 

-7^o($o(</'r+l,0r+2,A),$2(0r+l,0r+2, A),A) = 

by (3.47) and (3.65). 

To finish the proof, we proceed by induction. So, we assume that 

Jo(*i(0r+l,0r+2,A),$j+2(0r+l,0r+2,A),A) =0 (3.66) 

for < A; < j < r. Then 

Qk{^k{4>r+l,4>r+2,X),^k+2{(l)r+\,(l)r+2,Xj) =0 (3.67) 
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for < fc < j by theorem 3.16. Clearly, (3.50), (3.66), and (3.67) imply that 

P(4>o(</'fe+l, 0r+2, X))Tj+i{^j+i{(j)r+l,(pr+2, A), $j+2 (0r+l , 0r+2 , A) , A) = 0. 

(3.68) 

But, since 

^l + l (*J + l('/'r+l, 0r+2, A), $j+3((/)r+l, (/>r+2, A), A) G Y^'^^^^^^^^^^^x) 

by (3.63)-(3.64), we must, in fact, have 

-^i+l + 1 ('/'r+l, 0r+2, A), $j+3((/)r+l, (^r+2, A), A) =0 

by (3.47) and (3.68), and the proof is complete. □ 

In light of remark 3.14, the following corollary is a direct consequence of the 
above theorem. 

Corollary 3.19. Suppose s e Z>o, 3 < p < oo, and the maps , the neighborhoods 
K+i C X^'P C W^'P{B,R^) andNs+2 C C W^'P{B,M?), and Xq are as in 

proposition 3.16. Then for each {(j>s+i,(ps+2, ^} £ A^s+i x x (— Ao,Ao), the 

initial data 

{cl)\t=0,dt\t=0(l>) = i^o{ct>r+lAr+2,X),^l{^r+l,<f>r+2,X)) G 0'+''" X O'+^'P 

satisfy the compatibility conditions to order s. 

4. Local well-posedness 

In this section, wc prove local well-poscdness for the system (2.23a)-(2.23c) using 
the approach of [19]. The system under consideration is an initial-boundary value 
problem of elliptic-hyperbolic type, due to the presence of the equation (2.23b) 
in the system, and hence the results of [19] do not apply directly. However, the 
techniques of [19] are readily adapted to include non-local terms, and we will present 
an outline of the proof of this fact below. To conclude this section, we will apply 
the resulting existence theorem to establish the existence of dynamical solutions to 
(2.23a)-(2.23c). 

4.1. Setup and notation. For ease of reference, we adopt the index and other 
notational conventions of [19], with some exceptions, as pointed out below. We are 

interested only in the case n = 3, but it is convenient to treat the case of general 
n > 3. The number of components of the system of equations is in the case of 
elasticity equal to A'' = 3, but the treatment below applies to general N. 

Let < i, fc < n, 1 < a, ,5 < n, 1 < j, I < N. We work in a coordinate 
system (xi) and let t = xq. The summation convention is used. Further, we shall 
denote the unknown field by u rather than (f). Let Du = {dtU,dxiU, . . . ,dx„u), 
and DxU = [dx^u, . . . ,dx„u). Fix some s > n/2 + 1, s an integer, and consider 
the system in the domain f2 with boundary F of regularity class C^~^^, and denote 
Qt = [0, T)xQ and Ft = [0, T) x F where < T < oo. 

Following [19], we consider equations of the form 

dxiFj{t, X, u, Du) — Wj[t, X, u, Du] in fix, (4.1a) 

VjFj{t, X, u, Du) = gj(t, X, u, Du) in Ft, (4.1b) 

u = uq, dtu ~ ui in {0} x fl (4.1c) 

where, in contrast to [19], Wj[t, x, u, Du] is a functional of {u, Du) that we allow to 
be non-local. The properties required of w are specified below in Assumption w. 
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Assumption 1. We assume wo G W^^+H^^); «i ^ W''{^l). We assume -F, g G C+'-^U) 
where U is a. neighborhood of the graph of uo, ui, DxUq as in [19, §1], and 



We decompose hji = hjf + h^j^ where h^^f ~ ^{ji) ^^"^ symmetric part and /i"/^ = 
K^jT^ is the antisymmetric part. We assume that the symmetric part /i^^^ is of the 
form /i^j^;^ = O'^hji where 0{t,x,u, Du) is a vector field which is tangential to Ft 
and satisfies ^^'^ = 1, and /(,_,; is symmetric. 
For the elastic body, we have 

h^i = (4.2) 
while a'ji can be calculated in terms of the elasticity tensor L, cf. (2.16). 

Assumptions 2-5. The structural relations of the elastic body imply the hyperbolic- 
ity of the system. In particular, for the self-gravitating elastic body, the symmetry 
and coerciveness assumptions, Assumptions 2 and 3 of [19], hold. Further, Assump- 
tion 4 of [19] on the time components a'jf follow directly from the structure of the 
elastic system. For a discussion of the compatibility conditions on the initial data, 
see Assumption 5 in [19]. 

Following [19, p. 25], we introduce the spaces E"^, G^, and with norms 

\ 1/2 

\\dtu{t)\\i,s+,-A , 



and 




sup ||u||B-*(i), 

tl<t<t2 



Jtl 



respectively. 



Assum.ption w. For non-local w, we make the following further assumptions. We 

assume w to be well defined if the graph of (?i. Du) lies in a suitable subset of U, 
where U as in Assumption 1 above. There, the following conditions are imposed. 

(1) If 

u G no<j<,c^{[o,T],w'+'-^ift)), 

we have 

w[u,Du] € no<j<,C^{[0,T],W'-\^l)). 

(2) We require the map u h-)- w[u, Du] to be Lipschitz in the above topology. 

In particular, see section 4.3.1 below, we shall make use of the estimate 

\\d^{w{vi) - wivomL- < c||(at(«i - vo)\\e-. (4.3) 

(3) Finally, we require the following uniform estimate 

\\w\\e^ <cil + \\D^u\\L^)il + \\u\\Es+i) 

where c is a constant depending on HiJuHioo as well as the coercivity con- 
stants K, fi for the system. 
We have the following result, which is the analog of [19, theorem 1.1]. 
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Theorem 4.1. 

(1) Existence, regularity: There exists a unique < to < T, and a unique 
classical solution u E C^i^u, UTjJ of (4.1) with D''u{t) G L'^{Q) if < 
a < s + 1. Here D'^u denotes all derivatives of order a. 

(2) Continuous dependence on initial data. 

(3) Blow up: to is characterized by the two alternatives: either the graph of 
(u, Du) is not precompact in U or 

rt 

|D^u(T)||i:,oo(Q)dr — >■ 00, as t ^ to- 
la 

4.2. Linear systems. For a solution to tlie system (analogous to [19, (2.10)]) 

a^,(a'^9^,u) = w + d^J' in Ot, (4.4a) 

Va{a''''d,,u - D = h'''d,,u in Ft (4.4b) 

with the coercivity and structure conditions as in [19, §2], and in particular, the 
regularity assumptions (see [19, Assumption 2s, p. 26]) 

a,hGG' nC\nT) dta,dthGF\ (4.5a) 

w e G'-^ nW''^HlO,T],L^{n)) nW'-^'^{[0,T],W^{n)), (4.5b) 

f e G',Df' e W'^\[0,T],L^{n)), (4.5c) 

uoeW+^in) ui€W%n), (4.5d) 

we have the estimate 

||'"||B=+i(i2) < C^||u(ti)||£;3+i + ||w||g»-i 

rt2 



\\dMmLHn) + \\dmmLHa)dt (4.6) 
Jti / 

where 

c = c(k,^, \\a\\G''ncHnT)'\\9ta\\F-,\\h''\\G-nCHnT)'\\dth''\\F-)- 

Note that the system given in (4.4) is of a restricted form with 5 = 0, = 0. 
These terms can be absorbed into the others, cf. the discussion in [19, §2]. This 
is achieved by introducing the modified coefficients a^'',P,w as in [19, p. 31], see 
also (4.12) below. 

Applying the above estimate to a system of the form 

d^,{a'''d^,u) = w + d^J' in Qt (4.7a) 

Va{a"''d,,u - D = h>'d,,u + g in Tt, (4.7b) 

we get, in view of the above discussion, the inequality 



\\u\\Es^iit2) < d[ \\u{ti)\\Es+i + \\w\\gs-i + ||/||g= + ||5||g= 

''\\dMmL^n) + \\dt+'f{mLHn)dt 

1 

l-t2 

+ 



' Wtg\\w-(n) + \\dt+^g\\L^^n)ds^ (4i 
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4.3. Proof of theorem 4.1. In this section, we discuss the main steps in the proof 
of theorem 4.1. In the following calculations, we suppress the indices uj on u and 
the corresponding indices on aji,Wj, etc. First, we apply a time derivative to (4.1) 
which gives^ 

dx^a^'^dl^^tu) = dtw + d^J' in ^t, (4.9a) 

v^ia'^'^dl^^^u - n = h'^dl^, + g in Ft, (4.9b) 

where 

fit, X, u, Du) = a^^dl^^^u - dtF\ 9 = dt9- h'^dl^^.u. 

For the clastic system, the coefficients have no explicit time dependence, and the 
boundary condition is homogenous. So we have that 

.r=0, 5 = 0. (4.10) 

Next, the system is rewritten in the form 

dx,{a'''dl^^tu) =w + d^,r mnr, (4.11a) 

v^ia'^'^dl^.u - n = h^'^dl^.u in Ft, (4.11b) 

where 

f = f + v'g, a'^ = a''' - v'K^^ + w'^/i"' (4.12a) 
and 

w = dtw- {{d, + divw)/i"'=)a2^_,u + {da, v'^h:^')dl^^^u - d^g - {diYv)g. (4.12b) 

Remark 4.2. By introducing the modifications a^'',P,'W of a^'',P,'W as in (4.12), 
the resulting system (4.11) has no term g and also /i" = 0. Thus it is of the form 
of the system (4.4) considered in [19, theorem 2.4]. 

For the elastic system, we have f = g = h = and hence 

a''' = a''', w = dtW. (4.13) 

We have that a, / depend on it.x.u, Du) and il! depends linearly on D^u. 

For technical reasons, we assume s > n/2 + 2 for the differentiability index s. 
This is one more degree of smoothness than one would normally expect to require 
for a solution of a non-linear wave equation using energy estimates. However, this 
assumption reflects the use of f = dtu as the main variable in Koch's approach [19], 
which has one less degree of differentiability compared to u. The stated result for 
initial data with s > n/2 + 1 is recovered by a smoothing and a limit argument, see 
the discussion in [19, p. 33]. 

Next, we let Yt^h be the subset of 

Hr = ni<K,+iW^*'°°([o,T], w^^+i-^O)) 

of functions v that satisfy 9(^(0) = Ui, where Ui are the formal time derivatives of 
M for < i < s at t = 0, and < R. By choosing R sufficiently large, we can 

make sure this set is non-empty. 

The construction of solutions for the system (4.1) makes use of a standard fixed 
point argument, where one proves boundedness in a high norm and contraction in 



''We write d^^ ^ where dx^t 's used in [19] 
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a low norm. The high norm in this case is W^^^, which for the hnearized (time- 
diffcrcntiated) system corresponds to W^. The low norm for the time-differentiated 
system is W^. 

This type of argument has been carried out for a quasi-linear elliptic-hyperbolic 
system with no boundary conditions in [4] . The difference between that system and 
the present situation is that we have neumann-type boundary conditions, and the 
system has symbol depending on Du, i.e. it is fully non-linear. 

In the rest of this section, wc consider the details of the contraction estimate 
and the continuation property. The proof of the continuous dependence given in 
[19] can be readily adapted to the present case with the details given below. 

4.3.1. Contraction estimate. Define J G C{Yr^R, Ht) as the map v i-^ u, where u 
solves the linear system'' 

d^,{a'''{v)dl^^tu)=wiv) + d^,f in Or, (4.14a) 

Vc,{a'^''{v)<^,^^,u - riv)) = h\v)dl^^,u in Tt, (4.14b) 

dtu = ui, d^u = U2, in {0} x Cl (4.14c) 

where we have denoted a{v) = a{t, x, v, Dv) etc. 
For the case of the elastic system, we again have 

a''==a'*=, w = dtw, f = 0, h = 0, 

and in particular, 

w = DuW.dtu + Do^uW.Dxdtu 

where we have used D to denote the Frechet derivative. There is no explicit {t, x) 
dependence for this case. 
Next, we observe that 

dtw e G'-^, dt+^w G F° 

by Assumption w. From this, we see that for large R, the image of J lies in Y 
provided t is chosen small enough. 
We have 

FcC3([0,r],L2)nC72([0,r],W^i)n(:7i([0,T],W2-) 

One shows that the set Y is compact in the topology of the space defined by the 
right hand side of this expression. 

Let uq = J{vo), ui = J{vi). In order to derive the system solved by ui — uq, we 
let Z(u) = dxia^''{t,x,u,Du)dl.^ ^u, and note that Z{u\) — Z{uo) = ^Z{uo + 
X{ui — Uq)). Expanding this out gives 

Z(«i) - Z{uo) = S,.([ afdX\d',^^t{ui - uo) 
Jo 

+ dx,{l [Sl';„.(«i - UQ)]dl^^^{uQ + \{m - uo)))dX 
Jo 

+ Wlou-iDm - Duo)]52 ((mo + A(ui - uo)))d\ 

Jo 



typo in the boundary condition in [19, eq. (3.2)] is corrected here 
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where^ 



-ik 



0'\,Du 



a* (t, X, ua + \{ui — uo), Duq + \{Dui — Duq)), 
da. 
du 
dd 



if^u = ^(*'^''"o '^("i ~ uo),Duo + X{Dui - Duo)), 



dDu 



{t, X, Uo + X{ui - Uo), Duo + X{Dui - Duo)). 



Now, let 



::ik 



W 



ifdX 



f 



Jo 

iB{vi) — 'w{vo) 

- d^iif [df^u-{ui-uo)]dl^4uo + X{ui-uo)))dX 
Jo 

+dxA[ [atou-iDui - £'uo)]5^„tK + A(ui - uo)))dX 
Jo 

f{vi) - f{vo). 



Then ui — uo solves® 



vaia^'d'.^^tiui - uo) - n = -uo) + g 



(4.15a) 
(4.15b) 



where h, g can be calculated along the same lines as above. 

In particular, for the elastic system, we may calculate -w using w = dtw, h = 
f = g = 0. Note that g is non- vanishing in general due to contributions from d. 

We have the estimate"'^'^ 

\\dtw{t)\\L2 + \\dJ\\E^{t) + \\dtg\\E^{t) < c\\dt{vi - vo)\\e- (4.16) 

where we made use of 

w = w(v\) — w{yo) + terms involving a. 

As discussed above, w can be estimated given estimates for dtw, and hence dtw can 
be estimate in terms of df{w{vi) — w{vo)- The terms involving d can be estimated 
in terms of m — uq and hence can be absorbed when applying Gronwall. 
From Assumption w, (4.3), we have 

\\d'^{w{vi) - w{vo))\\l^ < c\\dt{vi - vo)\\e^- 

The required contraction estimate is obtained by applying the estimate for the linear 
system given in [19, theorem 2.4], as discussed in section 4.2, to the system (4.15). 
One checks that after suitable modifications, cf. section 4.2, the assumptions of 
[19, theorem 2.4] holds for this system, and this provides the needed contraction 
estimates. 



typo in [19, p. 32] is corrected here. 
®This corrects a typo in [19, p. 32]. 

^"^This corrects a typo in [19, p. 32]., which leads to an incorrect estimate for dt{ui — uq)- 
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Applying the inequality (4.8) with s = 1 to the system (4.15), we get the in- 
equality 



\\dt{ui - uo)\\e2 < c\^\\ui{0) - uo{0)\\e2 + \\w\\go + ||/||gi + IIsIIg 
+ £ \\dtw{a)\\L2 + liatfllsi + \\dt9{a)\\Erda^ 
use that ui,uo have the same initial data, and that w,f,g vanish at f = 

<c( [ \\dtwia)\\L2 + WdtfWE^ + \\dtg{a)\\E^da 



\J0 

use (4.16) 

< c\\dt{vi -vo)\\e'^ 
where we made use of the fact that w;, /, p all vanish at t = 0. 

4.3.2. Continuation principle. We next consider the proof of the continuation prin- 
ciple. Suppose the graph of {u, Du) lies in a compact set U. We need an estimate 
of the following form, cf. [19, eq. (3.3)] 



|w(t)|U.+i <^(u,J^ \\DMr)\\L^dT^ (l + lkol 



This estimate is proved by applying operators Dp of order s, tangent to Ft, to 
both sides of the equation (4.1), and applying the estimate for the linear system. 
Let s > n/2 + 1, and let u E 6"*+^ be a solution of the system (4.1). We then have 
u e r]Q<j<s+iCmo,TlW'+'^-^). One finds, cf. [19, p. 34], that Dj,u solves an 
equation of the form 

do,^{a'''d^kD-},u) =id + d^J' in fix, (4.18a) 

Viid^'d^.D^u - f) = h^'^d^^D-'pU in Tt, (4.18b) 

D'pu{Q) = u\, dtD'pU = u'^ (4.18c) 

where Mf,M| are obtained by formal calculations from the initial data uq,ui. Here 
a transformation which absorbs the terms g and /i"*^ has been applied, along the 
lines discussed in section 4.2, see [19, p. 34] for details. 

The basic energy estimate for systems of this type, see [19, theorem 2.2], gives 
an estimate of the form^^ 

\\u\\e.+^ < c [\\u\\es^.{Q) + ^* \\w\\l<t) + \\dtni^.{T)dT^ (4.19) 

As shown in [19], the L"^ norms in the right hand side of (4.19) that involve local 
expressions can be estimated at a fixed time in terms of 

(l+||I>2«||Loo)(l + ||u(f)||£. + l) 



ir 



The norms || • ||2 in [19, p. 35] should be || ■ ||^2. 
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with a constant depending on ||_Du||loo as weU as the coercivity constants k, /x. The 
term w contains DpW, and thus we need the nonlocal term w to satisfy at a fixed 
time an estimate of precisely this form, namely 

\\D'pw\\l. < c(l + \\D^u\\l^){1 + \\u\\es+i) 

which we can state as 

\\w\\e' < c(l + \\D\\\l'^){1 + \\u\\es+i). 

This estimate holds by Assumption w. 

4.4. Application to the elastic system. We are now ready to apply the local 
existence theorem 4.1 to our system (2.23a)-(2.23c). 

Theorem 4.3. Suppose s e Z>3 1 < p < oc, and let (^'|t=05 9t|t=o^') = {<pli,<p\) G 

Qs+2,p ^ Qs+i,p ^ vr*+2(g) X W+^iB) be the initial data from Corollary 3.19. 
Then there exists a to > and a unique classical solution (j)^ € C^{Btg U dBt^j) of 
(2.23a)-(2.23c) with D''(l)'{t) e L'^{B) for allO<\a\<s + l andO<t< to- 

Proof. First, we observe by Corollary 3.8 that non- local function A' ((^) = -6'^ ff-dAU{(f)) 
satisfies Assumption w of section 4.1 for (j/ in the open set C"*"^ C W^~^^{B). Since 
the initial data 

{<l>"\t=o,dt\t=o<l>") = i<l>h,<l>\) G d'+^'Pxd'+^-P C 0^+2.2x0^+1.2 ^ w''+^B)xW''+\B) 

from Corollary 3.19 satisfies the compatibility conditions to order s, the proof then 
follows directly from theorem 4.1. □ 
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Appendix A. Function Spaces 

A.l. W'^'P spaces. Give a finite dimensional vector space V, an open subset C 
with a C°° boundary. A; e Z, and 1 < < oo, we let W'''P{n,V) denote the 
standard Sobolev space for maps w : f2 — > V. If V = K, then we will just write 
W^'PiB), while if J? = 2 we set W''{n, V) = W^'^in, V). 
For these spaces, we recall the following results: 

(i) Differentiation 

" ^ ■ ^'''''(^) W'-^^Pifl) A =1,2,3 (A.l) 

defines a continuous linear map. 

(ii) (Multiplication Inequality) If 1 < p < oc. fci + fc2 > 0, ki,k2 > ks, fcs < 
ki + k2 — 3/p, then there exists a C > such that 

for ah u G W''^'P{n) and v G W''^'P{n). 

(iii) Letting B'''P{dQ) to denote the Besov spaces on the boundary dCl, the trace 
map 

W'''P{n) 9 M ^ tvonu G B''-^/P^P{dn) (A.3) 
is well defined and continuous for 1 < p < oo, and k — l/p> (cf. theorem 
7.70 in [1]). 
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A.2. Wg''^ spaces. For k e Z>o, 1 < p < oo, and (5 e M, we use Ws-''{R^,V) to 
denote the weighted Sobolev spaces for maps u : M'' — >■ V as defined in [5]. For 
open sets C R'^ for which \ 17 is bounded, we denote the restriction of the 
weighted spaces to these subsets by W^'^'^(ri, F). Fohowing [5], the negative index 
spaces '^(R^) for k G Z>o, 1 < p < oo, and 5 G R are defined by duahty. 
We rccaU the following facts about the weighted Sobolev spaces: 

(i) Differentiation 

dA : W^'^ift) Wtl'^i^) (A.4) 

defines a continuous linear map. 

(ii) (Sobolev Inequality) If k > 3/p, then there exists a C > such that 

hhfiny) < C|KllvKj='''(n,y) (^•^) 

for aU u e Wg-^iR^, V). Moreover, u e C?(R^y) and u{X) = o{\X\^) as^^ 
\X\ 00. (Se lemma A.7 in [27].) 

(iii) (Multiplication Inequality) If 1 < p < oo, ki + k2 > 0, ki,k2 > ^3, ks < 
ki + k2 — 3/p, and Si + 62 < Ss, then there exists a constant C > such 
that 

Il"^llw43,p(a) < C\\u\\^k,.,^^^\\v\\^,,.,^^^ (A.6) 

for ah u e Wg^^'^'in) and v G Wg^'^in). (See lemma A.8 in [27].) 

(iv) For 1 < p < 00, — 1 < (5 < 0, and fc e R, the Laplacian 

A = S'^^dAds : W^+^'^'iR'') W^:_p^{R^) (A.7) 
is a linear isomorphism with inverse given by the formula 

[A->(.)l(X) = -i-/^_^d»K (A^8) 

(The proof of this statement follows from using the fact that A : W^^^'^(M^) 
Wg^2i^^) is an isomorphism for k G Z>o and —1<S<0 (cf. [5, proposi- 
tion 2.2]) together with duality and interpolation^^.) 

For bounded O C R^, we let Eq denote an extension operator that satisfies 

Pa(w)||^^.P(M3,;.) <Cfc,p||w||^..p(o,^) VwG W'='f(f2,y), (A.9) 

and 

[d%Ea{u)] \n = d^u V u G W'''P{fl, V), \a\ < k. (A.IO) 



l^Herc \X\ = ^SabX'^XB. 

^"^As noted by Maxwell [25], the weighted spaces VV^'''(R"^) and their fractional extensions 
correspond to the spaces /ippj;,_^j_3 in [34, 35] (cf. remark 2 and theorem 2 in [35]). The 
following duality and interpolation results follow from remark 2, and theorems 2 and 3 in [35]: 

(a) For 1< p < 00 and 9 = p/{p - 1), l''^ (IR» ) is the dual of Wg''\M:^). 

(b) If 1 < pi < 00, 1 < p2 < 00 < e < 1, fc = (1 - e)ki + 9k2, 5 = (1 - e)5i + 95-2, and 
l/p = (l-6»)/pi +0/P2, then Wg'''(M.^) is the interpolation space [W^l'^iS.'^), W^^'''(M.'^)]e. 
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